Cancellation of nonlinear Zeeman shifts with light shifts 



K. Jensen/' 2 'Qv. M. Acosta,^ J. M. Higbie,"^ M. P. Ledbetter,^ S. M. Rochester,^ and D. Budker^' ^'0 

^ Niels Bohr Institute, University of Copenhagen, DK 2100, Denmark 
^ QUANTOP, Danish National Research Foundation Center for Quantum Optics 
^Department of Physics, University of California at Berkeley, Berkeley, California 94720-7300 
"^Department of Physics and Astronomy, Bucknell University, Lewisburg, Pennsylvania 17837 
^Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720 

Nonlinear Zeeman (NLZ) shifts arising from magnetic-field mixing of the two hyperfine ground- 
states in alkali atoms lead to splitting of magnetic-resonance lines. This is a major source of 
sensitivity degradation and the so-called "heading errors" of alkali-vapor atomic magnetometers 
operating in the geophysical field range {B ^ 0.2 — 0.7 G). Here, it is shown theoretically and 
experimentally that NLZ shifts can be effectively canceled by light shifts caused by a laser field of 
appropriate intensity, polarization and frequency, a technique that can be readily applied in practical 
situations. 

PACS numbers: PACS. 07.55.Ge, 32.60.+i, 42.65.-k 



I. INTRODUCTION 

Alkali-vapor atomic magnetometers [T] operating in 
the geophysical range of magnetic fields, where the Zee- 
man effect is very close to linear, are nevertheless quite 
sensitive to the nonlinear corrections to the Zeeman ef- 
fect (NLZ) which cause broadening and splitting of the 
magnetic-resonance lines, as well as line-shape asymme- 
tries that depend on the orientation of the sensor with 
respect to the field [H, @]. Thus, NLZ is responsible for 
sensitivity degradation and systematic (heading) errors 
in atomic magnetometers [1, [5| . Recently, collapse and 
revival of ground-state quantum beats associated with 
NLZ was studied theoretically [6] and experimentally Q, 
and a scheme for mitigating the effects of NLZ in an 
atomic magnetometer based on double-modulated syn- 
chronous optical pumping was realized Alternative 
approaches explored recently include the use of multi- 
quantum transitions and high-order atomic polarization 
moments 0, [H, [ol, [lo|- Here, we introduce an alterna- 
tive technique where NLZ shifts are compensated by light 
shifts due to an additional light field of appropriate in- 
tensity, polarization and frequency, which is possible due 
to the identical tensor structure of the splitting caused 
by the two effects. The present technique is free from 
some shortcomings of the alternative techniques such as 
complexity of implementation and/or degradation of the 
signal. Moreover, in contrast to the double-modulated 
synchronous pumping technique of Ref. the present 
approach works well when the ground-state polarization- 
decay rate approaches the NLZ frequency splitting (a 
common situation for practically important Rb and Cs 
magnetometers) . 

Compensation of NLZ with AC Stark shifts was re- 
cently investigated in the field of quantum information. 



theoretically as part of a quantum memory protocol [uj 
and experimentally for improving atomic spin-squeezing 
[l2[. The similarities of light shifts and Zeeman shifts 
was used in Refs. [l3|, where magnetic fields were 
simulated by light fields. It has also been proposed to 
use light shifts to measure parity violation in Fr atoms 
[l5[. In general, light shifts are important in precision 
measurements in atomic physics, for instance in atomic 
clock experiments, where one can use the famous "magic 
wavelengths" to eliminate light shifts, see for example 
Refs. [W, Il3 and references therein. 



II. THEORY 

For an alkali atom with nuclear spin /, total angular 
momentum F = / ± 1/2 and projection Mf of the to- 
tal angular momentum on the direction of the magnetic 
field, the Zeeman energies of the magnetic sublevels are 
to second order in the field B given by 
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where /i^ is the Bohr magneton, Ahfs is the hyperfine in- 
terval, and we neglect small corrections proportional to 
the nuclear magneton and the anomalous magnetic mo- 
ment of the electron. The first-order shift leads to preces- 
sion of the atomic polarization around the magnetic field 
with Larmor frequency vl = 2/iB^/[^(2/ + 1)], while 
the second-order shift leads to splitting of the magnetic 
resonances. For the F = 2 hyperfine manifold of ^^Rb, 
there are three AMp = 2 resonances with adjacent- 
frequency splittting 
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For ^''Rb atoms placed in a magnetic field of 0.5 G, we 
have = 350 kHz and S = 72 Hz. 
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Consider now ^^Rb atoms illuminated with linearly po- 
larized light near-resonant with the Dl transition. Based 
on the calculations in Refs. we can find analytical 

results for the AC Stark shifts of the energy levels in the 
F = 2 ground-state manifold. For light polarized along 
the magnetic field, we find a change in the splitting of 
adjacent AMj? = 2 resonance frequencies of 
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where eo is the permittivity of vacuum, Ap^^^ is the 
5^^1/2 excited-state hyperfine-structure coefficient mea- 
sured in Hz, Tdi is the natural linewidth of the 5^Pi/2 
excited state, A is the transition wavelength, Eq is the 
electric field amplitude of the light and the light fre- 
quency detuning A is measured relative to the F = 2 ^ 
F' = 1 transition and assumed to be much larger than 
the Doppler width and the upper-state hyperfine inter- 
val. For light polarized transverse to the magnetic field, 
we find that the change in splitting of the AMp = 2 reso- 
nances are of half the magnitude and of the opposite sign 
compared to the longitudinal case. In general, if the L2 
light is linearly polarized at an angle to the magnetic 
field, the differential shifts of the Mp sublevels are pro- 
portional to 3cos^^ — 1, see, for example. Problem 2.11 
in Ref. [lO]. Since the splitting increases in response to 
light of transverse polarization and decreases in response 
to longitudinal polarization, NLZ can only be compen- 
sated on this transition using longitudinal polarization. 
Using Eqs. ([2|) and ([3]) we can calculate the Dl light 
power needed to cancel the nonlinear Zeeman effect, and 
we find 
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Here we asssume that the atoms are kept in a cylindrical 
cell of diameter and that the average light intensity / = 
ceo£^o/2 ~ ^/ where c is the speed of light in 

vacuum, determines the amount of compensation. This is 
the case for atoms in an antirelaxation-coated cell, where 
each atom samples the light field in the entire cell volume. 
Similar calculations for light near-resonant with the D2 
transition shows that the change in splitting is of the 
opposite sign compared to the Dl case, and that roughly 
8-10 times more power is needed to cancel NLZ. 

In addition to shifting the magnetic resonances, light 
also broadens them. One source of this broadening is ab- 
sorption of light by the ground-state atoms. The broad- 
ening depends on the light power needed to compensate 
NLZ. For large detunings the broadening at the compen- 
sation power will be independent of detuning, since the 
power needed to cancel NLZ goes as A^ and the absorp- 
tion rate goes as 1/A^. Based on the calculations in 
Refs. [isl, [l^, we find the broadening at the compensa- 
tion power for Dl light with longitudinal polarization to 
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FIG. 1: (Color online) Schematic of the experimental appa- 
ratus. AOM - acousto-optical modulator; LP - linear polar- 
izers. A vapor cell with ^^Rb is located within a multi-layer 
magnetic shield. Coils within the inner shield is used to apply 
the Earth-range magnetic field and to compensate gradients. 
The intensity of light from laser LI is modulated as shown in 
the upper inset, so the same laser beam is used to optically 
pump and probe the atoms. Laser L2 is used to induce light 
shifts. 
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For D2 light with transverse polarization we find that 
the broadening at the compensation power is around 20 
times larger compared to the Dl case. It will therefore 
be preferable to use Dl light to cancel NLZ. 

Another source of light broadening results from the 
nonuniform intensity profile of the light beam over the 
cross-sectional area of the cell. If, as in our case, the 
compensating light illuminates only a small portion of 
the cell, the atoms are subject to short periods of com- 
pensating light at random intervals, so that the phase of 
each atom's evolution undergoes a random walk. This 
leads to dephasing of the atomic evolution at a rate 
Tdephase = ^'^ / 'T i whcrc (j) is the phase shift an atom 
experiences in each pulse of compensating light, and r 
is the average time between pulses. In order that the 
light properly compensates the nonlinear Zeeman effect, 
we must have (j) = 5t. The average time for an atom 
to cross the cell is on the order of d/v^ where v is the 
rms atomic velocity, and the probability that an atom 
passes through the light beam in one trip across the cell 
is on the order of h/d^ where h is the beam diameter, so 
T^d^/{hv). Thus 
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By inserting typical experimental values {B = 0.397 G, 
d ^ 2 cm and b ^ 1.5 mm), we find F^^ ^ 0.3 Hz and 
Fdephase ^ H Hz, and we see that dephasing dominates 
the broadening. 
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FIG. 2: (Color online) An example of data showing optical 
rotation of the probe light. At time = ms, the pumping stage 
ended. The signal undergoes several collapses and revivals 
during the decay. The magnetic field was set to 0.477 G. 
The inset shows a zoom-in revealing the fast oscillation at 
2uL = 667 kHz. 



FIG. 3: (Color online) Splitting between adjacent magnetic 
resonances (see Fig. U) as a function of L2 light power for 
polarizations along (circles) and transverse (squares) to the 
magnetic field (denoted tt and y polarization in the figure). 
Lines are linear fit to the data. The magnetic field was set to 
0.397 G. 



III. EXPERIMENTAL SETUP, PROCEDURE 
AND RESULTS 



Figure [T] shows the experimental setup. A parafiin- 
coated cylindrical glass cell (with diameter and length 
of 2 cm) at room temperature containing ^^Rb is located 
within a multi-layer magnetic shield with a system of coils 
in the inner volume designed to produce a homogeneous 
Earth-range magnetic field along the x direction. The 
atoms in the F = 2 ground-state are optically pumped 
and probed using a beam from a diode laser (LI) with 
frequency close to either the Dl or D2 line. The in- 
tensity of the light is modulated with an acousto-optical 
modulator (AOM). The same laser beam, initially lin- 
early polarized in the z direction, is used to create atomic 
alignment (the rank- two atomic polarization that is char- 
acterized by a preferred "alignment axis") and to probe 
the alignment via optical rotation due to the polarized 
atoms |2l|. A measurement consists of two consecutive 
stages. During the first, pumping stage, the laser in- 
tensity is square- wave modulated (duty cycle 1/8) at a 
rate equal to twice the Larmor frequency, and ground- 
state atomic alignment is synchronously pumped for 3 
ms. The beam diameter is ~ 3 mm, and the power dur- 
ing the "on" part of the cycle is ~ 5 mW. Once the atoms 
are pumped, the AOM is set to transmit around 4 /iW of 
light, which probes the atomic alignment created during 
the pumping stage. The optical rotation of the trans- 
mitted light is measured by a balanced polar imeter. An 
example of a free-induction decay (FID) signal (averaged 
over 1024 cycles) is shown in Fig. [2l The FID signal os- 
cillates at 667 kHz (see inset) corresponding to twice the 
Larmor frequency. The NLZ-induced beats in the signal 
are observed as a change in the oscillation amplitude with 
time. The signal has an offset of around 7 mrad due to 
imperfect balancing of the polarimeter. 

We model the measured FID signals with a function of 



the form 

S{t) = {C2 cos {27r[2iyL -S]t- a) ^Ci cos {AniyLt) 

+C2 cos (27r[2^L + S]t + a)} e'''^ , (7) 

which consists of a central component of amplitude C\ 
oscillating at twice the Larmor frequency and two side- 
bands with frequencies 2i^j^±(5 of equal amplitude C2 and 
opposite phase ±a. The three frequency components de- 
cay with time, and for simplicity it is assumed that the 
1/e decay time T is the same for all three components. 
We analyze the envelope of the signal by postprocessing 
the data with a digital lock-in amplifier with reference 
frequency 2vi,. The envelope R(i) can be calculated from 
Eq. ([7j) and we find 
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It is found experimentally that the demodulated FID sig- 
nals are well described by Eq. ([8]). For the data presented 
in Fig. [21 a fit of the envelope to Eq. ([8]) gives the value of 
5 = 65.37(7) Hz, consistent with the expected value due 
to the NLZ effect for the magnetic field oi B = 0.477 G 
used for the data presented in Fig. [2l The overall 1/e de- 
cay time extracted from the fit is 20.8(7) ms, limited by 
magnetic-field gradients and drifts in the bias magnetic 
field combined with signal averaging [for comparison, at 
lower magnetic field of 24 mG the relaxation time was 
measured to be 55(2) ms]. 

In order to cancel NLZ, linearly polarized light from a 
second diode laser (L2) tuned close to the Dl resonance 
was directed through the vapor cell along the z direction. 
The splitting of the AM^? = 2 resonances was measured 
as described above for different detunings, polarizations 
and powers of the L2 light. Figure [3] shows the splitting 
as a function of light power inside the cell for two polar- 
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FIG. 4: (Color online) Power spectrum of a FID signal show- 
ing the splitting of the AMp = 2 magnetic resonances due to 
NLZ (solid line) , and of a signal where NLZ was compensated 
by light shifts induced by the L2 light (dashed line). The 
magnetic field was set to 0.397 G. 

izations, along and transverse to the magnetic field (de- 
noted TT and y polarization in the figure), for a detuning 
A = -3.8 GHz from the F = 2 ^ = 1 Dl transition. 
The splittings were fit to straight lines yielding the slopes 
-9.4(13) Hz/mW and 5.6(8) Hz/mW, which within the 
uncertainties confirms that the light shifts due to y po- 
larized light have the opposite sign and are of half the 
magnitude compared to the light shifts due to tt polar- 
ized light. It should be noted that when the splitting is 
smaller than the linewidth of the resonances, it is difficult 
to extract the splitting from the FID signal. Therefore, 
in order to find the light power where the light shifts can- 
cel NLZ, an extrapolation from the linear fit was used. 
In Fig. [3] the light power where the splitting is zero is 
4.8(7) mW. 

Figure 2] (solid line) shows the power spectrum of a 
FID signal with the L2 light blocked. Three AMp = 2 
resonances are seen, split from each other in frequency 
due to NLZ. When the L2 light is on (dashed line), with 
the appropriate power needed to cancel NLZ, the peaks 
are combined into a single, stronger resonance. However, 
the L2 light also broadens the resonance. The FWHM 
of the resonances are 14.0(4) Hz (solid line) and 27.8(4) 
Hz (dashed line), giving a broadening of 13.8(6) Hz due 
to L2 light. By comparing the two spectra in Fig. [H 
it is seen that there is a shift of 8 Hz in the central 
frequency of the resonances. The shift is not thought to 
be due to the L2 light, but instead due to drifts in the 
current supply for the bias magnetic field in between the 
two measurements. 

The L2 light power needed to cancel the NLZ effect was 
measured for different frequencies of the L2 light close 
to the Dl resonance; the results are plotted in Fig. [5] 
(top). For these measurements, the magnetic field was 
set to B = 0.397 G, producing a NLZ shift S ^ 45 Hz. 
Also plotted in Fig. [5] (top) is the theoretical calculation 
for the compensation power given by Eq. (|4]), showing 
a reasonable agreement between theory and experiment. 
The broadening of the magnetic resonance due to the L2 
light at the compensation power was also measured and 
is plotted in Fig. [5] (bottom) together with the estimated 
broadening due to dephasing. In the frequency range 
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FIG. 5: (Color online) Top: L2 light power needed to cancel 
NLZ at a magnetic field of 0.397 G as a function of L2 laser fre- 
quency. Dots represent measured values, and the solid line is 
the theoretical calculation given by Eq. (|4|. Bottom: Broad- 
ening of the magnetic resonance due to the L2 light at the 
compensation power. Dashed line is the 11 Hz estimate due 
to dephasing given by Eq. (|6|). 

within -4 to 5 GHz of the F = 2 ^ = 1 Dl transition 
the broadening is approximately constant and consistent 
with the dephasing estimate. For the detuning A = 7.6 
GHz (which is on the F = 1 ^ F' = 2 resonance) a 
slightly larger broadening was measured. We note that 
the broadening due to dephasing can be avoided if the 
compensating light has a homegeneous intensity profile 
over the cell volume. In the experiment the beam diam- 
eter of the compensating light was rather small, and we 
therefore expect that it is possible to significantly reduce 
the broadening by for instance expanding the beam. 



IV. DISCUSSION AND CONCLUSIONS 

To be useful in practical magnetometers, our method 
for compensating NLZ should be robust against, for in- 
stance, intensity fluctuations of the L2 laser light. At 
Earth-range magnetic fields, NLZ splits the magnetic res- 
onances considered in the experiment by approximately 
70 Hz. We compensate NLZ by overlapping the reso- 
nances. The resonances will overlap if they are positioned 
within their width, which in our particular case is around 
30 Hz (see Fig. |4]). Since the splitting between the reso- 
nances is linear in L2 laser intensity, the relative intensity 
noise of the laser should be (much) less than the ratio of 
the linewidth to the splitting without compensation, in 
this case 30 Hz / 70 Hz 40 %. In practical situations 
one can easily stabilize lasers to have less than 1 % in- 
tensity noise. The intensity noise has therefore very little 
effect on the compensation. Similarly, we estimate that 
the frequency drifts of the L2 laser within easily achiev- 
able ±lMHz lead to sub-1-% changes in the induced light 
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shifts, much smaller than typical resonance width. 

In conclusion, using AC Stark shifts, we compensated 
the nonlinear Zeeman effect for ^^Rb atoms located in 
a magnetic field comparable to the Earth's magnetic 
field. The method can directly be applied to alkali- vapor 
atomic magnetometers in order to reduce heading errors 
and increase sensitivity. 



the manuscript and to E. Corsini for helpful discussions. 
This work has been supported by NURI grant HM1582- 
08-1-0006, ONR MURI and STTR grants. 



ACKNOWLEDGMENTS 



The authors are grateful to E. S. Polzik for encour- 
agement and support, to W. Gawlik for comments on 



[1] D. Budker and M. V. Romalis, Nature Physics 3, 227 
(2007). 

[2] V. Acosta, M. P. Ledbetter, S. M. Rochester, D. Bud- 
ker, D. F. Jackson Kimball, D. C. Hovde, W. Gawlik, 
S. Pustelny, J. Zachorowski, and V. V. Yashchuk, Phys. 
Rev. A 73, 053404 (2006). 

[3] S. J. Seltzer, P. J. Meares, and M. V. Romalis, Phys. 
Rev. A 75, 051407(R) (2007). 

[4] E. B. Alexandrov, Physica Scripta T105, 27 (2003). 

[5] Unpublished results by S. M. Rochester, J. Higbie and D. 
Budker show that heading errors are of less importance 
for alignment-based magnetometers compared to the tra- 
ditional orientation-based alkali- vapor magnetometers. 

[6] E. B. Alexandrov, M. Auzinsh, D. Budker, D. F. Kimball, 
S. M. Rochester, and V. V. Yashchuk, J. Opt. Soc. Am. 
B 22, 7 (2005). 

[7] E. B. Alexandrov, A. S. Pazgalev, and J. L. Rasson, Op- 
tika i Spektroskopiya 82, 14 (1997). 

[8] A. I. Okunevich, Optika i Spektroskopiya 91, 193 (2001). 

[9] V. M. Acosta, M. Auzinsh, W. Gawlik, P. Grisins, J. M. 
Higbie, D. F. J. Kimball, L. Krzemien, M. P. Ledbetter, 
S. Pustelny, S. M. Rochester, et al., Optics Express 16, 
11423 (2008). 

[10] V. V. Yashchuk, D. Budker, W. Gawlik, D. F. Kimball, 
Y. P. Malakyan, and S. M. Rochester, Phys. Rev. Lett. 



90, 253001 (2003). 
[11] T. Opatrny, Phys. Rev. A 74, 043809 (2006). 
[12] T. Fernholz, H. Krauter, K. Jensen, J. F. Sherson, A. S. 

S0rensen, and E. S. Polzik, Phys. Rev. Lett. 101, 073601 

(2008). 

[13] C. Y. Park, J. Y. Kim, J. M. Song, and D. Cho, Phys. 

Rev. A 65, 033410 (2002). 
[14] M. Vengalattore, J. M. Higbie, S. R. Leslie, J. Guzman, 

L. E. Sadler, and D. M. Stamper-Kurn, Phys. Rev. Lett. 

98, 200801 (2007). 
[15] M. A. Bouchiat, Phys. Rev. Lett. 100, 123003 (2008). 
[16] Y. Jun, H. J. Kimble, and H. Katori, Science 320, 1734 

(2008). 

[17] V. V. Flambaum, V. A. Dzuba, and A. Derevianko, Phys. 

Rev. Lett. 101, 220801 (2008). 
[18] W. Happer and B. S. Mathur, Phys. Rev. 163, 12 (1967). 
[19] B. S. Mathur, M. Tang, and W. Happer, Phys. Rev. 171, 

11 (1968). 

[20] D. Budker, D. Kimball, and D. DeMille, Atomic physics 
An exploration through problems and solutions (Oxford 
University Press, 2008), 2nd ed. 

[21] W. Gawlik, L. Krzemien, S. Pustelny, D. Sangla, J. Za- 
chorowski, M. Graf, A. O. Sushkov, and D. Budker, Appl. 
Phys. Lett. 88, 131108 (2006). 



